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ABSTRACT 

We present simulations of dry-merger encounters between pairs of elliptical galaxies with dark matter 
halos. The aim of these simulations is to study the intergalactic stellar populations produced in both 
parabolic and hyperbolic encounters. We model progenitor galaxies with total-to-luminous mass ratios 
Mt/Ml = 3 and 11. The initial mass of the colliding galaxies are chosen so that M1/M2 = 1 and 10. 
The model galaxies are populated by particles representing stars, as in Stanghellini et al. (2006), and 
dark matter. Merger remnants resulting from these encounters display a population of unbounded 
particles, both dark and luminous. The number of particles becoming unbounded depends on orbital 
configuration, with hyperbolic encounters producing a larger luminous intracluster population than 
parabolic encounters. Furthermore, in simulations with identical orbital parameters, a lower Mt/Ml 
of the colliding galaxies produces a larger fraction of unbounded luminous particles. For each modeled 
collision, the fraction of unbounded to initial stellar mass is the same in all mass-bins considered, 
similarly to what we found previously by modeling encounters of galaxies without dark halos. The 
fraction of intergalactic to total luminosity resulting from our simulations is ~ 4% and ~ 6% for 
dark-to-bright mass rations of 10 and 2 respectively. These unbounded-to-total luminous fractions are 
down from 17% that we had previously found in the case of no dark halos. Our results are in broad 
agreement with intergalactic light observed in groups of galaxies, while the results of our previous 
models without dark halos better encompass observed intracluster populations. We suggest a possible 
formation scenario of intergalactic stars. 

Subject headings: Galaxies: elliptical and lenticular, cD; interactions; stellar content. Stars: AGB and 
post-AGB. Planetary nebulae: general. 
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1. INTRODUCTION 

The population of gravitationally unbounded inter- 
galactic (IG) stars has been recently observed in dif- 
ferent environments, such as in poor groups (Feldmeier 
et al. 2004, Castro- Rodriguez et al. 2003), compact tidal 
groups (White et al. 2003), and nearby galaxy clusters 
(Ferguson et al. 1998; Durrell et al. 2002; Feldmeier 
et al. 2004, Gerhard et al. 2005), as well as in higher 
redshift clusters (Zibetti et al. 2005, Zibetti 2008). The 
amount of observed IG starlight seems to correlate with 
the environment, with more IG starlight found in pop- 
ulated clusters, and toward the cluster centers. Zibetti 
(2008) studied 683 clusters of galaxies, and found that 
the fraction of diffuse to total starlight decreases with 
the distances of the unbounded stars from the galactic 
centers, a result that is also supported by Murante et al. 
(2004) 's predictions. 

Several mechanisms have been proposed to explain the 
origin of the IG starlight. Muccione & Ciotti (2004) as- 
sociated the presence of IG stars to the evaporation of 
galaxies due to tidal interactions with other cluster mem- 
bers, and with the cluster potential well. Conroy, Wech- 
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sler & Kravtsov (2007) showed that more than 80% of 
the stars from tidally disrupted sub-halos contributes to 
the IG starlight. 

As an alternative to tidally produced intracluster stars, 
we (Stanghellini et al. 2006, hereafter Paper I) have 
shown that dry mergers could produce an IG stellar pop- 
ulation compatible with what has been observed to date 
in the intracluster. Our models of Paper I show that 
during the merging stage a number of stars from the 
progenitor systems become unbounded, and the IG to 
initial starlight obtained was broadly consistent with the 
observed intracluster star counts. 

Murante et al. (2007) analyzed a large number of clus- 
ters from a cosmological simulation. They find that most 
of the stars in the IG stellar population come from the 
formation of the brightest cluster galaxy, and that most 
stars of this diffuse component come from merging events 
and only a minor fraction of the stars derive from tidal 
interactions. 

In this paper we extend our models of interacting galac- 
tic pairs, whose components have dark matter halos. 
There is now ample evidence that disk galaxies should 
host a fair amount of dark matter to account for their 
kinematics (see, van Albada et al. 1985 and references 
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thereafter). A further indication that elliptical galax- 
ies have dark halos comes from cosmological simulations, 
and the hierarchical merging paradigm (see, e.g. Onorbe 
et al. 2007). It is thus appropriate to repeat the simu- 
lations of colliding pairs of elliptical galaxies as in Paper 
I, this time assuming that the initial galaxy models have 
dark halos. We perform our simulations keeping all other 
assumptions unvaried with respect to those of Paper I, 
so we can directly compare the results from both sets of 
simulations. 

In Section 2 we introduce our models and all the as- 
sumptions made for the galaxies, the dark halos, the stel- 
lar populations, and the collisions. Section 3 illustrates 
the results from the dry-merger simulations, with partic- 
ular attention to the production of the IG stellar pop- 
ulation. The discussion, including comparing the new, 
dark-haloed models with those of Paper I, is in §4. Sec- 
tion 5 presents the conclusions of our study. 

2. MODELS 

We have performed a series of simulations of dry merg- 
ers of pairs of spheroidal galaxies with dark matter halos. 
We chose the mass ratio between the two colliding galax- 
ies to be M 2 /M x =l and 10 1 . 

2.1. Galaxies, Initial conditions 

The initial conditions of the galaxies forming the col- 
liding pairs are similar to those described by Gonzalez- 
Garcfa & van Albada (2005), who performed a study 
of the encounters between two spherical systems with 
dark matter halos. We use the isotropic spherical Jaffe's 
(1983) approach to model the initial conditions of the lu- 
minous matter, and the Hernquist (1990) model for the 
dark matter halos. 

The potential for a Jaffe model is: 
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where G is Newton's constant, and rj is the half-mass 
radius of the luminous component. 
The potential of the Hernquist model is: 
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where a is the scale length used. The half-mass radius 
of the dark halo component is equal to (1 + y/2)a. We 
combined the Jaffe and the Hernquist model to reproduce 
the haloed galaxy. We obtained a two-parameter family 
of galaxies that depends on the luminous-to-dark mass 
ratio, Ml/Mh, and on the half-mass radius ratio, 



rL 1/2 /r Hl/ 



(l + \/2>' 

The distribution function of the combined system is 



(3) 



h(E)=h(E) + f H (E). 



(4) 



1 Hereafter the suffix 1 refers to the less massive galaxy of 
the merging pair, l refers to the luminous mass, h refers to the 
dark halo mass and x to the total mass of an individual model, 
Mj = Mjj + Ml- Unless otherwise noted, M is the total mass. 



An algorithm yielding fh(E) and fn(E) along the two- 
parameter family was developed by P.A.H. Smulders & 
M. Balcells (1995, unpublished, see Gonzalez-Garci'a & 
van Albada 2005 for a detailed description of the algo- 
rithm). 

In the following, we adopt non-dimensional units with 
G = 1 for Newton's constant of gravity. The theoretical 
half-mass radius of the Jaffe model, rj, and the total 
luminous mass of the fiducial galaxy in each run are also 
set to 1. The models may be compared with real galaxies 
using the following scaling, similar to the one employed 
in Paper I: 



[M] = 4 x 10 11 M 0) 
[L] = rj = 5 Kpc, 
[T] = 8.34 x 10 6 yr. 
By adopting these units, the velocity unit is: 

\v] = 587 km s -1 . 



(5) 
(0 
(7) 
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The initial systems are spherical and non-rotating, 
with an isotropic velocity distribution. In this regard 
they are similar to the set of spherical isotropic systems 
studied in Paper I. 

The observed amount of dark matter in elliptical galax- 
ies is controversial. Elliptical galaxies lack the reliable 
kinematic tracers that allow the detailed dynamical stud- 
ies of spiral galaxies. Recently there have been improve- 
ments in recovering the presence of dark matter in ellip- 
tical galaxies, such as in the X-ray study by Humphrey et 
al. (2006), via gravitational lensing (Treu & Koopmans 
2004, Mandelbaum et al. 2006), and through stellar dy- 
namics (Thomas et al. 2007). Chakrabarty & Raychaud- 
hury (2008) have applied globular cluster kinematics to 
estimate a ratio of ~10 for the dark to luminous mass 
in the elliptical galaxy NGC 4636. Romanovsky et al. 
(2003), from observations of the kinematics of planetary 
nebulae in the outer parts of elliptical galaxies set a firm 
upper limit to the amount of their dark matter, although 
it has been noted that their results may be model depen- 
dent (Dekel et al. 2005). 

In the present models, we have initially used a range 
of Mh/Ml values, and calculated the circular velocity 
curve for those realizations (bearing in mind that we use 
a composite of the Jaffe and Hernquist models), aiming 
at a flat-topped curve for a significant radial range. We 
obtained the best match for Mh/A/l = 2 and 10, and a 
halo scale length a — 2 which yield component ratios of 
Mt/Ml=3 and 11. These ratios are consistent with the 
observations of Mandelbaum et al. (2006). The models 
presented here can be viewed as extrapolations of the 
Paper I models, where no dark matter halo was included 
(i.e., Mt/Ml=1). Dark halos have been populated with 
5x 10 5 particles for all galaxy models, except the smallest 
galaxy in run lOhP which has ten times less particles. 

The projected surface density of luminous material of 
such models presents a slope that decreases roughly as 
i? 1 / 4 , which makes it a suitable representation for ellip- 
tical galaxies, although the central parts present a cusp. 

In Table [1] we summarize the characteristics of the ini- 
tial models. Column (1) gives the run identification code, 
where the initial number indicates the mass ratio of the 
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TABLE 1 

Galaxy input parameters 



Run 


M T /M L 


M2/M1 


n 


«i 


b 


E or b 


lhP 


11 


1 


40 


0.316 








lhZ 


11 


1 


40 


1.048 





0.250 


lOhP 


11 


10 


19.48 


0.336 








lhP3 


3 


1 


40 


0.316 









interacting galaxy pair, the small letter denotes the im- 
pact parameter (h is for head-on impact), and the cap- 
ital letter denotes the energy of the orbit, with P for 
parabolic and Z for zero energy at infinity (or hyper- 
bolic) orbit; column (2) gives the total-to-luminous mass 
ratio of the initial galaxies; column (3) gives the mass 
ratio of the colliding galaxies (where the masses repre- 
sent the total masses for each galaxy); column (4) gives 
the initial separation r^, in model units; column (5) gives 
the initial relative velocity of the galaxy pair, Uj, also in 
model units; and columns (6) and (7) give respectively 
the impact parameter, and the orbital energy of the ini- 
tial setup. 

To appropriately scale the models with mass other than 
1, we impose that the luminous mass of initial systems 
must lie on the fundamental plane of galaxies which, ac- 
cording to J0rgensen et al. (1996), gives: 

logL ~ 0.78 log M T + constant. (9) 

and we use Fish's (1964) law to scale the radii of the 
luminous and dark components of the two systems. 

Model lhP is an equal mass encounter (Mi = M2) with 
the interacting galaxies set on a parabolic orbit. The cen- 
ters of the two galaxies are placed at an initial distance 
of 4 Ri, where R% = 10 rj is the galaxy radius. Model 
lhZ is an equal mass encounter with the progenitors on a 
mildly hyperbolic orbit, with all the other parameters as 
in model lhP. Model lOhP is an encounter between two 
galaxies with a mass ratio of 10 on a radial parabolic en- 
counter, initially placed at a distance 3 Ri + R2, where 
Ri and R2 are respectively the radii of the less and the 
more massive galaxies. Following the definitions above, 
R 2 = 10 rj and Ri = 3.16 rj. Finally, Model lhP3 has 
the same orbital parameters as model lhP, but with dif- 
ferent total to luminous components ratio of the progen- 
itor system, Mt/Ml=3. 

The relative range of the mass of the galaxy pairs cho- 
sen for the simulations, Ma/Mi=l or 10, are the extremes 
of the observed merging galaxy sample by van Dokkum 
(2005). 

2.2. Stellar components 

The luminous components of our galaxies models have 
been populated by particle stars that follow Salpeter's 
(1955) initial mass function (IMF), ^f(rn) cx to~ 2 - 35 , and 
are distributed into three representative mass bins. A 
detailed discussion of this implementation is given in Pa- 
per I. The population of each mass bin corresponds to the 
integration of the Salpeter's mass function in that bin, 
scaled to the entire population considered, and ignoring 
stars outside that mass range. The luminous particles in 
our fiducial galaxy model have mass distributions such 
as that, after accounting for the Salpeter's IMF, the total 
galaxy luminous mass is equal to unity. 

In Table [2] we summarize the luminous population of 



TABLE 2 
Stellar input 

bin [Mq] <S> bin TSS parUcles 

[Mel 

0.85 - 1.4 0.5151 156060 3.3 x lO" 6 1.32 x 10 b 

1.4 - 3.0 0.3443 52121 6.6 x 10~ 6 2.64 x 10 5 

3.0 - 8.0 0.1406 7050 2.0 x 10~ 5 8.00 x 10 5 

total 1 215231 1 

the input fiducial galaxy. Column (1) gives the mass bin 
in solar masses, column (2) gives the mass fraction in 
that mass bin, column (3) gives the number of luminous 
particles in the bin, and column (4) the (dimensionless) 
mass of particles in each bin. 

In runs lhP, lhZ and IhPZ both galaxies have num- 
ber of particles per mass bin as in Table [2] The stellar 
content of Table [2] also represent the more massive galaxy 
in run lOhP, while the less massive galaxy (Ml = 0.1) 
has an initial set up with ten times less particles in each 
mass bin, but with the same mass per particle. 



2.3. Integration method and stability 

We have used the parallel treecode GADGET-2 
(Springel 2005) on the MareNostrum Super Computer 
at the BSC. Softening was set to ~ 1/10 of the half-mass 
radius of the smallest galaxy (e = 0.075). The tolerance 
parameter was set to 6 = 0.8. Quadrupole terms were 
included in the force calculation and the maximum time 
step was set to 1/50 of the half-mass crossing time. 

We have checked the stability of our input initial model 
for 33 half-light crossing times (see formula 2.40 in Bin- 
ney & Tremaine 2008, page 64; - 0.5 Gyrs). The test 
shows that the system relaxes for about 6 time units 
and remains stable thereafter. The inner parts present a 
shallower profile at the end of the relaxation time. This 
initial relaxation is due to the presence of the particle 
softening in the code. The stability runs are compa- 
rable to that performed for ten times less particles by 
Gonzalez-Garcia & van Albada (2005) with small differ- 
ences in the inner parts due to the smaller softening used 
in the present simulations. 

The merged system was evolved for more than 10 dy- 
namical half-light crossing times of the merged body af- 
ter merging, to allow the system to relax (reach virial- 
ization) . Conservation of energy is good in all the runs; 
energy variations are less than 0.5 per cent. 

3. RESULTS 

The evolution of the runs and the final states are quite 
similar to those already described by Gonzalez-Garcia & 
van Albada (2005), we refer to the interested reader to 
that paper for a detailed description of the merger stages. 
We present a brief description here. 

Galaxies in simulation lhP evolve towards a first pas- 
sage through pericenter at t — 80 (in model units of 
time). At this first encounter the systems exchange en- 
ergy, the halos of the two galaxies already merge, and a 
fraction of both dark and luminous matter gains enough 
energy to escape. A similar situation is found in simula- 
tion lOhP, although with different timescales. Galaxies 
in simulation lhZ ', which has the highest initial orbital 
energy, merge after the first encounter and thus there 
is only a single release of particles. Simulation lhP3 
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shows a different behavior due to the lower relative con- 
tent of dark matter. The two systems meet at t ~ 80, 
but the halos do not merge at this first encounter. A 
number of particles of the two systems, both dark and 
luminous, gain enough energy to escape. The two system 
meet again at t ~ 1 50 and finally merge, again exporting 
some particles (dark and luminous) to the IG popula- 
tion. The simulations are stopped at time t = 200, i.e. 
t ~ f .67Gyrs. 

Figure fT] shows the final distribution of the various 
types of particles in our encounter simulations; in each 
panel the x-axis represents the distance from the merger 
center, R, and the y-axis the binding, Ebmd, in dimen- 
sionless model units. We plot the results of different 
particles: from the left to the right panels: dark matter 
particles, then low-mass, intermediate-mass, and high- 
mass luminous particles are plotted. The dark and lumi- 
nous particles show similar qualitative behavior. When 
the galaxy systems pass through the pericenter, parti- 
cles gain kinetic energy and may end up having a pos- 
itive binding energy. Simulation lhZ presents a larger 
number of particles with positive binding energy than all 
other types of encounters. Simulations lhP and l/iP3 
show a moderate number of expelled particles. Finally, 
simulation lOhP shows very small number of particles 
with positive binding energy. 

In Figure [2] we plot the logarithm surface density cal- 
culated in annuli of increasing radius from the merger's 
centers versus R 1 / 4 , where R is the distance from the 
merger center. The different lines represent the total 
surface density (solid) and the surface density for each of 
the luminous populations (see caption). 

The intermediate regions present a profile close to a de 
Vaucoleurs profile (see full straight line). The innermost 
(i? 1 / 4 < 0.7) and outer regions deviate clearly from this 
profile. The inner parts are affected by the fact that our 
initial models have a core-like shape after the relaxation 
run. The outer regions deviate from the de Vaucouleurs 
profile due to the escaping IG material in a similar fash- 
ion as it was explained in Paper I. However, this devia- 
tion is less pronounced in the present models given the 
lower fraction of escaping particles. Deviations from the 
de Vaucouleurs' slope affected by the IG particles in the 
outer regions are important outside i? 1 / 4 ~ 2—2.5. These 
values are in agreement with those found in Paper I. The 
surface density excess with respect to a de Vaucouleurs 
fit (solid thin line) is larger for model lOhP than for 
model lhP at large radii. This is due to the deposition 
of most luminous particles of the small galaxy at these 
large radii. In all models and at all radii the different 
populations contribute similarly to the IG population. 

In each of the proposed simulations, a fraction of par- 
ticles escapes the main galaxian potential, in Table 3 we 
give the resulting IG population produced in each run; 
column (I) gives the run code; column (2) gives the mass 
bin; column (3) gives the fraction of unbounded luminous 
mass (with respect to total luminous mass); column (4) 
gives the fraction of the luminosity of unbounded par- 
ticles, L u , with respect to total luminosity Lt and (5) 
gives the fraction of the luminosity of unbounded parti- 
cles versus the luminosity of the merger remnant, L m . 
We also give, for each run and in column (3), the mass 
of the total (luminous and dark) unbounded mass as a 



fraction of the total initial mass. 

By examining Table 3 we note that the fraction of grav- 
itationally unbounded mass is higher in the total than 
in the luminous matter. This is occurring because the 
dark halos extend further out than the luminous parts of 
the galaxies, thus the presence of the dark halo poten- 
tial places the luminous matter deeper into the potential 
well. The fraction of total unbounded mass is roughly 
the same in runs lhP and lOhP, thus it is independent 
on the initial mass ratio. It is on the other hand higher in 
IhPZ due to the different internal structures and relative 
energies while varying the total mass of the dark halo. 
The situation is different for lhZ, where the larger avail- 
able energy results in 12% of the mass being expelled by 
the merging process. 

The luminous mass expelled from the encounters of 
galaxies with dark halos in Table 3 is less than what 
was found in the models without dark halos (i. e., Paper 
I): we find that the modeled stellar-to-total unbounded 
mass after the parabolic encounters of galaxies with dark 
halos is about 15 times smaller than in models without 
dark halos (it goes down to 10 times for the hyperbolic 
encounters) . 

Gonzalez-Garcia & van Albada (2003) and Nipoti et al. 
(2003) demonstrated that a few merger events between 
systems on the FP do not destroy such relation (although 
see Nipoti, Treu & Bolton 2009). The merger remnants 
in our sample must obey a law like log L ~ 0.78 log Mt 
+ constant (J0rgensen et al. 1996; equation 9). In the 
present study we have that Mt = Ml + Mjj. Now, if we 
introduce this in Equation 9, we get: 



1 Ll 78 1 ML1 

log — = 0.78 log-— 

L/2 Ml2 



0.78 log 



(M T i/M 
(M T2 /M 



L1J 



L2 



(10) 



We must translate the mass of our stellar particles into 
light. Our simulations are gravity-only simulations and 
we do not have their luminosity. To convert our stellar 
masses in luminosity we use the relation M T /M L cx M T 17 
from Hyde & Bernardi (2009) combined with equa- 
tion 9. Both the progenitor and the merged systems 
are on the FP relation, but the unbounded material is 
not. So, at a first step we want to compute the following: 
L u /L m = (Li + L 2 - L m )/L m = (Li/L m ) + (L 2 /L m ) - 1. 
Using the above mentioned formulas we get: 



L u _ (Ml_\0.94 , , M 2 NO. 94 

L m y mJ k mJ 



(id 



And we can derive the fractions, which arc included 
in Table [3l Finally, we find that the ratio of luminous 
unbounded to total material L u /Lt = L U /(L U + L m ) = 
1/[(L U + L m )/L u ) = 1/(1 + (L m /L u )), also included in 
Tabled 



4. DISCUSSION 

In reviewing all our merging models we find that he 
total amount of the IG starlight is determined mainly 
by two factors: the amount of dark matter in the col- 
liding galaxies, and the energy of the encounter. The 
first factor is clearly accounted for when we look at sim- 
ulations lhP and lhP3 from this paper and model lhP 



Intergalactic Stellar Population from Mergers of E's 



5 



r 



lhP 

Halo 



lhP 

Popl 



1 « 

500 1000 500 1000 ' 

R R 




500 1000 1 
R 



lhP 

Pop3 



500 1000 
R 




lhZ 

Halo 1 ™ 



200 400 

R 




400 




1 






1 






lhZ 




Pop3 







200 400 
R 




lOhP 

Halo 



100 200 

R 











o 




1 




H 


lOhP 


M 


Popl 


CO 













o 










o 














lOhP 


H 


lOhP 


N 


Pop2 


M 


Pop3 


co 




co 


1 




200 



200 



Fig. 1. — Binding energy vs. distance from the merger's center of models lhP (top row), lhZ (second row from the top), lOhP (third 
row from the top), lhP3 (bottom row). Columns contains, respectively from left to right, dark matter particles, low-mass population of 
luminous particles, intermediate-mass population of luminous particles, and high-mass populations of luminous particles. Units are in the 
internal system of units, see text for details. 
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Fig. 2. — Surface density vs. i? 1 / 4 for the different simulations. 
From top to bottom, respectively, we show lhP, lhZ, lOhP, and 
lhP3. The different lines are the results for the different stellar 
populations. Low-mass population: Solid line; intermediate-mass 
population: Dashed line; and high-mass population: dottcd-dashed 
line. A de Vaucoleurs profiles is shown as a thin solid straight line. 
Units are in the internal system of units. For details, see text. 



TABLE 3 

Results 



Run 


bin [M Q ] 


M u (%) 


L U /L T 


Lu/Lm 


(1) 


(M e ) 

(2) 


(%) 
(3) 


(%) 
(4) 


(%) 

(5) 


lhP 


0.85 - 1.4 
1.4 - 3.0 
3.0 - 8.0 


0.33 
0.34 
0.31 








total lum 


0.33 


4.37 


4.57 




Total 


2.06 






lhZ 


0.85 - 1.4 
1.4 - 3.0 
3.0 - 8.0 


2.07 
2.04 
1.87 








total lum 


2.04 


5.92 


6.29 




Total 


12.15 






lOhP 


0.85 - 1.4 
1.4 - 3.0 
3.0 - 8.0 


0.63 
0.60 
0.52 








total lum 


0.60 


2.38 


2.43 




Total 


2.41 






IhPZ 


0.85 - 1.4 
1.4 - 3.0 
3.0 - 8.0 


2.36 
2.22 
2.38 








total lum 


2.32 


6.17 


6.58 




Total 


3.25 







in Paper I : The fractions of IG to total luminosity are 
~ 4%, ~ 6% and 17% for dark mass amounts of, respec- 
tively, 10, 2, and times the amount of the luminous 
mass. On the other hand, increasing the mass ratio of 
the colliding galaxies (model lOhP) does not produce a 
significant increased IG population: The unbounded lu- 
minosity produced in this model is small when compared 
with the other simulations. This is reflecting the fact that 
most of the IG material originates, in this case, from the 
low-mass galaxy. This also implies that satellite-like or 
minor mergers should be at least between two and ten 
times more frequent to get a similar IG stellar population 
than that produced by major mergers. 

We can compare the present results and those in Paper 
I with other models and simulations. The semi-analytic 
predictions of Purccll, Bullock & Zentner (2007) claim 



that the intrahalo light should rise from 0.5 to 20% from 
galaxy halos to poor groups, while it should be between 
20 — 30% for rich groups and clusters. Such numbers for 
the IG light in poor groups arc in rough agreement with 
our results for the simulations presented here for models 
with a dark halo, while their numbers for clusters are in 
closer agreement with our results from Paper I. 

Observations might indicate an increase of the rela- 
tive IG to total starlight as the mass of the group in- 
creases. Loose groups show fractions of the order of less 
than 2% (Castro-Rodriguez et al. 2003) while, clusters 
like Virgo present fractions around 5-10% (Arnaboldi et 
al. 2003) and massive clusters fractions larger than 10% 
(Feldmeier et al. 2004), although Zibetti et al. (2005) 
find that the fraction of IG to total starlight is almost 
independent on the cluster richness. 

Aguerri & Gonzalez-Garcia (2009) show that a galaxy 
entering a dense environment, such as a cluster, it is 
stripped away from the outer parts (mostly dark mat- 
ter), due to tidal stripping with the potential well of 
the cluster as well as due to galaxy harassment. The 
harassed galaxy has, as a consequence, a lower dark- to- 
luminous mass ratio. This transformation is more se- 
vere in dense environments where faster encounters are 
likely. Thus, denser environments should host systems 
with lower dark-to-luminous ratios. Conversely, groups 
should host systems with larger dark-to-luminous ratios. 

In summary, in loose groups the potential well of the 
group is not large enough to strip the galaxy from all 
its dark matter halo. This scenario is likely to present 
mergers of galaxies with a fair amount of dark matter 
involved, like the ones investigated in this paper. The 
shielding effect of the dark halos should then prevent a 
large fraction of IG stars in these environments. Dense 
groups and rich clusters, on the contrary, are environ- 
ments where tidal stripping and tidal harassment with 
the halo potential should affect primarily the outer parts 
of the galaxies, stripping away the halo; this might have 
the effect of lowering the initial relative velocities and 
the components dark-to-luminous ratio, and the galaxy 
pairs might undergo encounters such as those presented 
of lhP3 in this paper, or those presented by Paper I. This 
could explain why these environments present a larger IG 
light fraction. 

The parameter space explored here is naturally limited. 
In this sense, this paper and Paper I must be viewed as 
a first attempt. We test the assumption that most of 
the IG population come from dry-mergers, so gaseous 
mergers involving disk galaxies, star formation, and re- 
processing are excluded from our modeling. Mergers be- 
tween gas-free disk (SO-like) galaxies have recently been 
explored (Wu & Jiang 2009) with results compatible with 
those presented here. 

Other caveats of our approach could be that we do 
not model the collisions arising from cosmological sim- 
ulations, and thus we do not have fully cosmological 
conditions. Our dark matter halo is a Hcrnquist pro- 
file, which is a cored halo profile. However, cosmologi- 
cal simulations tend to produce cuspier density profiles 
(see e.g., Navarro, Frcnk & White 1996). We have run 
a test simulation with an NFW dark matter halo pro- 
file, constructed with the software described by Widrow 
& Dubinski (2005), with initial conditions identical to 
run lhP. The fractions of unbound light are: L u /L m = 
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4.43% and L U /L T = 4.24%. These figures are roughly 
similar to those found in runs with shallower inner pro- 
files. Altough the NFW is deeper at the central parts, 
we must note that the unbound matterial is extracted 
from the outer layers of both colliding systems where the 
potential well is shallower. 

Also, in clusters the matter distribution in satellite ha- 
los is modified and would be better tested with other 
profiles like those proposed by Springel et al. (2008). 
Such problems should be explored in the future using 
a re-simulation of a cluster or group area from a fully 
cosmological simulation. Finally, we must note that the 
IG light could have been pre-processed in groups and 
later accreted to a cluster, contributing with their share 
of already unbound stars (see, e.g. Rudick, Mihos & 
McBride, 2006). 

5. CONCLUSIONS 

We have modeled dry-mergers between pairs of ellip- 
tical galaxies populated with luminous particles, sur- 
rounded by dark matter halos and predicted the ex- 
tent of the intergalactic stellar populations produced in 
these encounters. These simulations modeled different 
orbital parameters, mass ratios, and varying initial dark- 
to-luminous mass ratios. These simulations are an exten- 
sion of those presented in Paper I, which were without 
dark matter halo. 

The merger process injects orbital energy into the fi- 
nal system. Some particles, both luminous and dark, 
gain enough energy to become unbounded, populating 
the modeled intergalactic space. The resulting fraction 
of unbounded to initial dark matter is larger than the 



corresponding fraction of luminous particles. We also 
find that the presence of a dark matter halo results in a 
lower amount of free-floating luminous particles, trans- 
lating into an expectation of more intracluster or intra- 
group stars in low-dark matter galaxy encounters. 

An interesting point would be to know when particles 
become unbound and which dynamical process (violent 
relaxation, tidal stripping at pericenter...) is the domi- 
nant one. However, this is out of the scope of the present 
article and will be investigated in forthcoming papers. 
Finally, the fraction of IG starlight observed in the dif- 
ferent simulations presented in this paper are in rough 
agreement with those reported in analytical models, and 
in observations of poor galaxy groups. The fractions pre- 
sented in Paper I are in closer agreement with those ob- 
served in the intracluster. 

We conclude that dry merging is a viable process to ex- 
plain the IG unbound stellar population in many specific 
cases of low group density. 
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